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Ab stra ct 

The dynamics of flexible spacecraft are not usually well known before launch. This makes it 
important to develop controllers for such systems that can never be destabilized by perturbations 
in the structural model. Virtual passive controllers, or active vibration absorbers, possess this 
guaranteed stability property; they mimic a fictitious flexible structure attached to the true 
physical one. 

This report analyzes the properties of such controllers, and shows that their disturbance 
absorption behavior can be naturally described in terms of a set of virtual zeros that they 
introduce into the closed-loop dynamics of the system. Based on this analysis, techniques are 
then derived for selecting the active vibration absorber internal parameters, i.e. the gain matrices 
of such controllers, so as to achieve specified control objectives. 

Finally, the effects on closed-loop stability of small delays in the feedback loop are investigated. 
Such delays would typically be introduced by a digital implementation of an active vibration 
absorber. It is shown that these delays only affect the real parts of the eigenvalues of a lightly- 
damped structure. Furthermore, it is only the high-frequency modes that are destabilized by 
delays; low-frequency modes are actually made more heavily damped. Eigenvalue perturbation 
methods are used to obtain accurate predictions of the critical delay at which a given system will 
become unstable; these methods also determine which mode is critical. 
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1. Introduction 

An important practical difficulty in the control of flexible space structures (FSS) is that their 
dynamics are usually only poorly known. This is a consequence of the great difficulties involved 
in determining on-orbit behavior from pre-flight ground vibration tests. The sensitivity of an 
FSS control system to variations in the dynamics of the structure it is applied to is therefore a key 
practical consideration. A case that is particularly to be avoided is that of an unduly sensitive 
controller which gives good closed-loop performance for the nominal structure but destabilizes 
the true, off-nominal, vehicle. 

Such considerations motivate the important class of dissipative controllers [1]. These control 
systems, based on the use of compatible (physically collocated and coaxial) actuators and rate 
sensors, guarantee closed-loop stability regardless of uncertainties present in the structural 
model. A related type of controller that has recently been proposed [2][3][16][17] is that of 
virtual passive controllers, in which the control system mimics a fictitious (or virtual) flexible 
structure attached to the physical structure to be controlled. The controller design problem then 
reduces to that of choosing the natural frequencies, damping ratios and mode shapes of this 
virtual structure in such a way as to satisfy the specifications placed on the closed-loop system. 
Such compensators are ideally suited to a high-authority/low-authority control scheme; they 
provide the stabilizing low-authority inner loop, so allowing a higher-performance feedback 
technique with no guaranteed stability, for instance a linear quadratic regulator (LQR), to be 
applied safely. 

The first subject of this report is to develop methods for designing virtual passive controllers 
when the primary control objective is to achieve vibration absorption at various specified 
disturbance frequencies; the controller in this case is referred to as an active vibration absorber 
( AVA ). This control objective amounts to requiring that transmission zeros [5] at these 
frequencies be introduced into the transfer functions between disturbance inputs and measured 
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outputs. A necessary preliminary step is therefore an analysis of the properties of the zeros of 
AVA-controlled flexible structures. The results obtained will be shown to be natural extensions 
of the open-loop FSS zeros relations of [4]; they also help to quantify the additional flexibility 
offered by AVA control over that provided by state feedback. Furthermore, this zeros analysis 
illuminates the significant differences between the case where the AVA sensor/actuator pairs are 
coincident with the disturbance sources or primary output sensors, and that where they are at 
distinct locations. 

It should be noted that virtual passive controllers may be based on the use of compatible 
acceleration, rather than rate, measurements; this is the approach that will be taken here. The 
ease with which accelerometers can be distributed throughout a structure makes implementation 
of such a scheme extremely practical. Methods will then be derived for designing an 
accelerometer- based AVA which places the new virtual zeros at desired locations. The 
remaining design freedom is next quantified and exploited to satisfy some specified secondary 
control objectives, for instance minimization of the sensitivity of the closed-loop poles. 

The effects on closed-loop stability of small delays in the feedback loop are also investigated in 
this report. Such delays would typically be introduced by any digital implementation of an 
active vibration absorber, so analysis of their effects is of great practical importance. First-order 
generalized eigenvalue perturbation methods are used to study these questions, and certain 
fundamental results obtained. As a first point, it is shown that delays acting on a lightly-damped 
flexible structure will basically only affect the real parts of the closed-loop eigenvalues. 
Furthermore, not all modes are destabilized by delays. Only the high-frequency modes are 
destabilized; the low-frequency modes are actually made more heavily damped. Accurate 
predictions are obtained of the critical delay at which a given system will become unstable; these 
methods also show which mode is the critical one. It should be noted that this will not 
necessarily simply be the highest-frequency mode of the structure. In fact, numerical simulations 
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suggest that this is very rarely the case. Finally, all the results obtained in this report are 
illustrated by application to simple examples. Further information on the work reported here is 
also given in [22], 


2. AVA Equations 

Consider an n-mode model for the structural dynamics of a non-gyroscopic, non-circulatory FSS 
acted on by low-authority control inputs u, disturbances w and high-authority control/reference 
signals r. This model can be written as 

Afq + Cq + ATq = flu + fl 0 w + B h t, 

( 2 . 1 ) 

y = H a q + H r q + H d q, 

where the mass, stiffness and damping matrices satisfy M = M T > 0, K = K T > 0 and 
C=C T > 0, respectively. The generalized coordinate vector q is (n x 1), u is (m x 1), and the 
measured output vector y is (p x 1). 

A virtual passive controller, or active vibration absorber (AVA), for this system is a low- 
authority controller which takes the form of a fictitious flexible structure attached to the physical 
structure. This can be written in modal form as 

% + diag(2£ vi CQ vi )Ti v + diag(co 2 jT\ v = 3>Ju v , (2.2) 

where is (n v x 1) with n v * n in general. The input to this virtual system consists of m 
accelerometer measurements from the physical structure, collocated with the low-authority 
actuator stations. 


u v = B T q. 


(2.3) 


Accelerometers are small, accurate and inexpensive, making them extremely suitable for 
practical use. This is the reason for the selection of acceleration feedback here. It should be 
noted, though, that displacement feedback would have to be added if rigid-body modes were to 
be controlled also. The analysis required in that case is entirely analogous to that which follows, 
and so will not be given explicitly. 

The outputs from the virtual structure are m acceleration measurements, by analogy with the 
physical measurements, and collocated with its input stations, 

y v = <Mv (2 ’ 4) 

The manner in which the resulting low-authority control signals are applied to the physical 
system, i.e. the interconnection between the virtual and physical structures, is given as 

u = y v -Gu v , ( 2 - 5 ) 

where the positive semi-definite symmetric gain matrix G — is chosen so as to ensure 
stability of the closed-loop system regardless of perturbations in the FSS open-loop dynamics. 


Combining (2.1) and (2.2) gives the composite equations of motion 


r M 0\ 
0 I 


q + 


fc o 

v 0 diag(2£ vi co vi ) j 


L 

f K 0 > 

r 

v 0 diagicol )j 


■( 


Bu + B d w + B h 
<t> T u . 

v “v / 


( 2 . 6 ) 


where q = (q T i^) r . Using (2.3) - (2.5), the right-hand side can be rewritten as 
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(2.7) 


r B[y v -Gu v ]+B D y/+B H r y 

V j 

_ ( M> v ii v - BGB T q + B D w + B H 

" , ; 

Taking all terms involving x to the left-hand side then yields the closed-loop equations of motion 
Mij + C<} + jfq = B d w + B h t, 

(2.8) 

y = 77j + 77^ + /7,q. 


where 


M = 


[M+B® v <Z> t v B t ] -B® 


T dT 


c = 


(C 


-®[B 


0 


and 


0 diag(2C vl co vl ) 


(2.9) 


- (K 
K = 


0 '' 
2 


0 diag(0)J 


Similarly, B D ={B t d 0 f , H a =(H a 0), etc. To simplify the notation in what follows, we 
define the (n x n v ) matrix X = fiO v . The composite mass matrix then becomes 


M = 


[M + XX T ] -X 
-X T I 


( 2 . 10 ) 


Now, the composite stiffness and damping matrices are positive semi-definite, by inspection. 
Defining the gain matrix G as in (2.5) ensures that M is also actually positive definite, so 
guaranteeing that the composite equations of motion (2.8) represent the dynamics of a flexible 
structure. Positive definiteness can be proved as follows: for any q = (q q v ) we have 
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( 2 . 11 ) 


q r Afq = q T (M + XX T )q - 2q r Xq v + r\ T v l\ v 
= q r Mq + (X r q-n v ) r (X r q-il v ) 

= q r A/q + z T z > 0, 

as M> 0. Thus, the closed-loop poles produced by the AVA will be stable, regardless of any 
perturbations in the open-loop dynamics (2.1). It is interesting to note that the vector z is closely 
related to the control signal applied to the physical structure by the AVA: this is just, by (2.3) - 
(2.5), 

u = -O v [X r q-fi v ] = -<I> v z. (2-12) 

The AVA design problem to be studied in this paper can be stated as follows: select the 
parameters {<w vi }, {£„,}, <b v and B (i.e. the AVA sensor/actuator locations) so as to minimize the 
frequency response magnitude of the closed-loop measured outputs {>”;} at some specified 

frequencies. These frequencies will typically be those at which the harmonic disturbance inputs 
{w ; ) are expected to occur. This vibration isolation problem will be considered to be the primary 

control objective; any additional design freedom available will then be used to satisfy other 
secondary objectives, for instance minimization of the sensitivity of the closed-loop poles to 
perturbations in the plant. 

The vibration isolation problem is equivalent to requiring that transmission zeros [5] be 
introduced, at the specified disturbance frequencies, into the transfer functions of the composite 
system between specified disturbances and specified outputs. It will be shown below that the 
solution of this question depends centrally on whether the AVA sensor/actuator pairs can be 
coincident with either the disturbance sources or the measured output (y) sensors. In order to 
fully investigate the solution methods and properties in the two cases, they will be considered 
separately in the two sections that follow. 
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3. Coincident Case 

The properties of the zeros of a flexible structure controlled by an AV A have not been studied to 
date. This analysis will now be given as a necessary basis for the results to follow. It will also 
serve to illuminate the differences between the case where the AVA is coincident with the 
disturbances or measured outputs and the general non-coincident arrangement. 




Taking the Laplace transform of (2.8) gives the polynomial matrix representation [8] 

P(s)q(.y) = B d w(s) + B H r(s), 
y(s) = H(s)q(s), 


for the composite system, where P(s) = s 2 M + sC + K and H(s) = s 2 H a +sH r + H d . The 
transmission zeros of this system between disturbance input j and measured output i are those 
values of s for which the rank of the system matrix [5] 






-w 0 


(3.2) 


is reduced. This polynomial matrix definition is more convenient to deal with than the rational 
one involving the transfer function T {j (s) = b i {s)P(sy l b j \ so long as the input-output loop 

examined is completely controllable and observable, the results obtained are the same in both 
cases. 


These zeros can be studied by means of a generalization of the method employed in [4] [6] [9] to 
investigate the zeros of open-loop flexible structures. This makes use of the QR decompositions 
[7] of b, and hf, i.e. =(q M Q bl ){p b 0 0) r and hf = Q h r h 

= (q M Qhi){Ph 0 "■ Of, with Q b and Q h both orthogonal matrices and p b , p h non-zero 
scalars. Applying these transformations to S(s) allows it to be reduced to a canonical form from 
which the transmission zeros can be obtained by inspection. In particular, 
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(3.3) 


i n, 


As p b and p h are non-zero, the first row and column of this matrix may be omitted when 
searching for values of s which make it singular. Thus, the zeros are just those values for which 

(& ^ 
lo h 


S(s) = 


M Q o ? 


(3.4) 


has less than full rank. 


So far, collocation of AVA and disturbance has not been assumed. The effect of making this 

restriction is to obtain a transformed system matrix which is block triangular, leading to a simple 
decomposition result for the resulting zeros. In detail, if B = b ; then X = b ; <E> v , so Ql 2 X = 0. 

Thus, by (2. 10) we have 


(f "l 57 


ra, o', (qi 2 mq„ o' 


U 12 

0 


-XQ 


h2 


fez, owa, 

0^ 

and 


0 > 

to 

o 

o 

Cl 

o 

h 


lo 

h 

10 I) 


also, 

therefore indeed block lower triangular, of the form 

* (Q T b2 [s 2 M + sC+K]Q hl 0 A 
l -s 2 XQ h2 P v (s) )’ 


(3.5) 


are block diagonal by inspection. S(s ) is 


(3.6) 
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where P v (s) = s 2 I + s ■ diag{2£ vl CO vi } + diag{co The top left block of this equation can be seen 
to be just the polynomial matrix which defines the structural zeros [4] of the open-loop physical 
flexible structure, while the lower right block is the denominator matrix of the AVA virtual 
structure. 

The zeros of the composite system between disturbance j and output i are therefore made up of 
three distinct sets: 

(a) the sensor zeros, consisting of a single zero at the origin for a rate measurement, and two 
for acceleration; 

(b) the structural zeros of the SISO open-loop physical system for this disturbance/output pair; 

(c) the poles of the virtual system (2.2). These will be referred to as the virtual zeros of the 
composite system. 

The first two sets of zeros are independent of the choice of AVA parameters; they are just those 
zeros which the open-loop physical structure initially possessed. Adding closed-loop zeros at 
desired frequencies is therefore achieved by placing the poles of the virtual structure at these 
locations. This makes sense physically, as the virtual system will then resonate at these 
frequencies, so absorbing energy and minimizing vibration. It should be noted that this result is 
entirely independent of the parameters of the open-loop physical structure; it holds whenever the 
disturbance source and AVA sensor/actuator pair are coincident. It will be shown in the next 
section that the results that hold in the non-coincident case are not so simple; in particular, they 
dfl depend on the open-loop model. 


11 


It is well known [5] that the zeros of a given system cannot be altered in any way by state 
feedback. Thus, AVA control provides an additional degree of design freedom over that 
provided by state feedback. The virtual zeros are, in fact, closely related to the disturbance zeros 
defined in [10] for the case of dynamic output feedback. However, that paper did not examine 
flexible structure applications or properties. 

It was assumed above that the AVA was coincident with the disturbance source of interest. An 
entirely similar result applies if the AVA is instead coincident with the output sensor. (The 
transformed system matrix S(s) is then block upper triangular rather than lower.) The poles of 

the virtual system again become the additional zeros of the composite system, and the only 
significant difference appears when the problem is generalized to the multi-disturbance/multi- 
output case. An AVA which is coincident with disturbance j can be used to introduce a zero at 
this disturbance frequency for arbitrary output locations. Thus, such an absorber provides 
"global" attenuation of this particular disturbance. By contrast, an AVA coincident with output 
sensor i can only absorb the disturbance at this physical location; the other outputs will still be 
affected by it. However, a compensation for this "local" attenuation is the fact that the AVA can 
now absorb multiple disturbances at this output station. In particular, an AVA of order n v can 
isolate the output from ny disturbances, generally acting at different frequencies and applied at 
different locations. This is a very powerful property. 

We have therefore shown that, in the coincident case, the natural frequencies {co vi } of the virtual 
system should be chosen to be equal to the disturbance frequencies which are to be absorbed. 
Also, the absorber influence matrix B is constrained by the requirement that the AVA be 
coincident with either the disturbances or the outputs. Therefore, the only design variables 
remaining to be assigned are the damping ratios {£„,} and modal matrix <1> V . These may be 

selected as follows. 
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The level of virtual damping is chosen based on a trade-off which can be viewed in either the 
time or frequency domain. Decreasing the virtual damping gives a deeper trough in the 
frequency response of the composite system at the desired attenuation frequency; however, the 
trough is narrower, giving less robust vibration absorption in the event of perturbations in the 
disturbance frequency. Similarly, lower damping gives a smaller steady-state time response to a 

sinusoidal disturbance at the specified frequency, but at the expense of a longer decay time for 
the initial transients. An optimum choice for {£ VI } will therefore depend on the details of the 

particular application considered; as a general rule though, simulations to date indicate that 
values in the range 0.05 - 0.10 give good results. 

The virtual influence matrix d> v can actually be chosen arbitrarily without altering the virtual 
zeros at all. It can therefore be regarded as making up a set of design variables which may be 
used to satisfy some secondary system requirement, over and above the primary one of vibration 
isolation. Of course, for the case of a single degree-of-freedom AVA, altering d> v amounts to a 
simple scaling of the virtual output y v ; it can be shown to be dynamically equivalent to scaling 
the mass of the virtual system (2.2). In general, increasing O v increases the gain matrix G 
(equation (2.5)), and so the control effort u (equation (2.12)); it is therefore essentially a control 
gain-type variable. 

Two particular secondary design requirements will be investigated here: minimization of an 
'LQR-like' steady-state error quantity, and minimization of the sensitivity of the poles of the 
closed-loop system. In the first case, the objective function is of the form 

J =\\yssi+p\Kt (3 - 7) 

where the system is excited by the harmonic disturbance signal. The control weighting p is used 
as in the standard LQR method to allow the relative importance of control effort and output 








performance to be varied. In the second problem studied, the quantities to be minimized are the 
sensitivities of the poles of the composite system, i.e. the eigenvalues of 



0 / 


(3.8) 


The sensitivity of the eigenvalue A. of A is measured by its condition number c,, which is 
defined [11] as follows: if A, has corresponding right eigenvector x, (Ax^AjX,) and left 
eigenvector y f (y?A = A ( y"), with ||x ( .|| 2 = |y,|| 2 = 1, then its condition number is 



(3.9) 


This quantity, which can be seen to be always greater than or equal to unity, measures the 
magnification between a perturbation in A and the resulting induced perturbation in A,. 
Choosing <I> V so as to reduce the {c*} yields a closed-loop system which is robust as well as 

guaranteed stable; this is certainly desirable in practice, given the uncertainties inherent in FSS 
models. Examples of the use of both the methods just outlined for selecting <I> V will be given 

later in the paper. 



4. Non-Coincident Case 

It may not always be possible in practice to place the AVA sensor/actuator pairs at either the 
disturbance sources or the primary output sensors. The closed-loop properties that then apply 
are, as has already been noted, considerably more complicated than those that held in the 
coincident case. In particular, the transmission zeros introduced by a non-coincident AVA are no 
longer simply equal to its own poles; a degree of dependency on the physical structure enters into 
the problem of choosing the AVA natural frequencies. This model-dependency is not unduly 
serious, however, as closed-loop stability is still guaranteed regardless of open-loop 
perturbations. (See [1] for a similar observation concerning positive real controllers.) 


14 



't i 

w 


The analysis required in the non-coincident AVA case proceeds initially in the same manner as 
that given in the last section. Equations (3.2) - (3.4) apply entirely without change, and the extra 

complexity only arises when we apply the required transformations to the composite mass matrix 
M. As B * b; here we no longer have Q[ 2 X = 0, so M and S(s ) are not block lower triangular. 

Instead, if we define the quantities 


Y = Qlx 


and 


Z = Q T h2 X, 


(4.1) 


(2.10) implies that the transformed mass matrix becomes 


[Q T b2 MQ h2 + YZ T ] -Y 
-Z T I 



r Ql 2 0^ 

m[ 

Qh2 

0) 

f 


lo i) 

' 

< 0 

h 

\ 


(4.2) 


Consequently, 

S(s) = 


[QLp(s)Qh2+s 7 yz t ] - s 2 Y ' 
s 2 Z T P v (s) 


(4.3) 


where P(s) = s l M + sC + K . Now, the transmission zeros of the composite system are those 
values of s which make S(s ) singular. Thus, making use of the standard determinantal identity 
[ 12 ] 


< 3 - 


det{A} • det{D - CA B) 


(4.4) 


for partitioned matrices, we have that any zero s must either satisfy 

det(a r 2 P(j)a : +j ! KZ r ) = 0 


(4.5) 


or 
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det{P v (J) - s 4 Z T [Ql 2 P(s)Q h2 + s 2 YZ T y l Y) = 0. 


(4.6) 


Note that, in the coincident case, either Y or Z is zero; (4.5) then reduces to the defining 
equation for the structural zeros of the physical structure, and (4.6) to det{P v (.?)} = 0, the defining 

equation for the virtual zeros. The results of the last section are therefore indeed a special case of 
those obtained here. The fact that the structural zeros depend on the applied feedback in the non- 
coincident case is analogous to the observation [13][14] that the zeros of a flexible structure can 
be shifted by dampers only if these are not coincident with the primary sensors and actuators. 


The AVA design problem in the general case thus reduces to that of finding a P v (s) = s' I + 
s diag[2£ vi co vi } + diagiQ)l} which satisfies (4.6) for some desired closed-loop zero s. This 
equation becomes much simpler to solve when the AVA has a single degree-of- freedom ( n v = 1 ); 
Y and Z are then column vectors and P v (s) is a scalar, so (4.6) reduces to 


p v (s ) = j V[# 2 p( 5)0« + * VrV 


(4.7) 


All quantities on the right-hand side are known, while the left-hand side is just s 2 + 2£ v co v s + 0) 2 V . 
Thus, we have a relation of the form 


2£ v o)J + 6)1 = a + j(3. 


(4.8) 


where s = a+ jt. The required virtual frequency and damping ratio are then given from the 
known quantities a, (5, o and t as 


az-Ba 


and 


(4.9) 


c,= 


i 

2cos 
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These solutions always exist; however, they may be imaginary in certain cases, making the 
resulting P v (s ) physically meaningless. Numerical results to date indicate that such infeasible 

solutions occur only when high gains (i.e. large <I> V ) are used and s is positioned in the vicinity 
of a closed-loop structural zero. This choice of zero causes the matrix in (4.5) to be nearly 
singular, leading to numerical ill-conditioning problems when performing the inversion in (4.7). 
However, this case is very unlikely to arise in practice: if a structural zero is near the disturbance 
frequency, it will already provide sufficient attenuation. There will therefore be no need to place 
a new virtual zero there as well. 

The multi-degree-of-freedom AVA problem must be solved using an entirely different approach, 
which will now be briefly described. For simplicity, we shall neglect the small amounts of 

damping present in both the original flexible structure and the desired virtual zeros 
[s t : l = 1 The second matrix in (4.6), denoted by W t say, is then clearly real. The solution 

P v (s ) consequently proves to be undamped also, with natural frequencies which must satisfy 
det{diag(o) 2 vi ) + [diag(sf) + W t ]} = 0, (4.10) 

l = 1 These simultaneous determinantal equations may be readily solved using a symbolic 

manipulation package such as Mathematica [15], giving the required virtual natural frequencies. 
Of course, just as in the single degree-of-freedom case, these values may on occasion turn out to 
be imaginary, and so infeasible. The likelihood of this occurring in practice is a subject requiring 
further research. 

5. Closed-Loop Equations with Delay 

Introducing a delay into the feedback loop used to implement the AVA controller corrupts the 
symmetry of the composite mass, stiffness and/or damping matrices. This in turn prevents the 
unconditional, model-independent stability guarantees that held in the undelayed case. 
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In detail, if we assume a small delay X in the feedback signal, we can use the first-order 
exponential approximation 


e~ sr ~\-sT 


to compute the perturbed matrices of the closed-loop system. These are 


M t - M + 1 M, 


C T ~C + rC 


K r ~K. 


Note that the composite stiffness matrix is, to first order, unchanged by the introduction of the 
delay. The perturbation matrices are given as 

_ f 0 -BQ> v diag{2£ vi a) vi j\ .... 


0 -B<& v diag(CQ 2 J 


The important point to note is that neither of these matrices are symmetric. Thus, the closed- 
loop system obtained for a non-zero delay does noi correspond to a composite fictitious flexible 
structure. It therefore does not possess guaranteed stability, as will now be demonstrated in more 
detail. 

The delayed system has poles given by the roots of the characteristic equation 


det{A 2 A/ T + XC t + A"} = 0. 
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This equation is not amenable to solution in its present form, as it is in terms of a polynomial, or 
lambda , matrix. To obtain useful numerical or analytical results, it must be transformed to a 
simpler first-order form. Two possible main approaches exist: to reduce it to a generalized 
eigenvalue problem of the form 

Ax = Afix, (5.8) 


or to reduce it even further to a standard eigenvalue problem by inverting the matrix B. The 
generalized eigenvalue formulation is the option to be preferred here, as performing a matrix 
inversion corrupts and confuses the matrix perturbation expressions (3.2) - (3.6) that form the 
basis of the semi-closed- form results that will follow. Defining the state vector 



(5.9) 


yields a generalized eigenvalue problem (5.8) for the delayed system (5.7), where 
( -C x -K^ 


A = 


0 


and 


B = \ 


'M r <n 


X 

0 / 


(5.10) 


(5.11) 


In terms of the perturbation expressions given previously, these matrices can be written as 


A = A 0 + 5A 


and 

B — Bq + SB, 


where 

A) =1 


f' ? 

-r 

& 

II 

( -C 0^ 

l ' 

0 , 


,0 oj 


(5.12) 

(5.13) 

(5.14) 
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(5.15) 


B 0 = 


f M 0 ^ 
0 / 


and 5B-x- 


M 0 
0 0 


A recent perturbation result for the generalized eigenvalue problem (where such results are 
considerably more difficult to obtain than for the standard eigenvalue problem) can now be 
applied to these matrices. Suppose that x is the right eigenvector of the nominal generalized 
eigenproblem corresponding to eigenvalue X, i.e. the solution of 


A)X = XB 0 x, 


(5.16) 


and y is the corresponding left eigenvector, i.e. the solution of 

y"A = y"XB 0 . (5-17) 

Then a first order approximation to the eigenvalue of the perturbed generalized eigenproblem 
(5.8), with A and B as given by (5. 12) and (5.13), is [18] 

X + 8X =y H (A) + 8A)x / y H (B 0 + 5B)x. (5.18) 

Now, premultiplying (5.8) by y H yields 

X = y H A qX / y H B 0 x. (5.19) 


Thus, we can rewrite and expand (5.18) as 

(A + 8X)y H (B Q + 5B)x ~ Xy H B 0 x + Xy H 8Bx + <5Ay w fl 0 x 

- y w (A> + <5A)x 

= Ay w fl 0 x + y h 8Ax. (5.20) 

The perturbation in the generalized eigenvalue X can therefore be written, to first order, as 

8X~[y H 8Ax-Xy H 8Bx]ly H B 0 x. (5-21) 

This result is valid for any sufficiently small matrix perturbations. For the present application, 
we can exploit the special form of Eqs. (5.14) and (5.15) to further simplify it. As a first step. 
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since the matrices SA and SB are both proportional to the delay x, so the eigenvalue perturbation 
will be also. This clearly makes sense. Furthermore, we can write the eigenvectors x and y of 
the nominal system in block form as 


x = 


and 


y = 


y„ 

jb. 


where these sub- vectors satisfy 


x fl = Ax*, 

[X 2 M + XC + K]x b = 0 

-y*K=Xyl 
yf[A 2 M + AC + tf] = 0. 


( 5 . 22 ) 


(5.23) 


(5.24) 


(5.25) 


Eq. (5.22) then reduces to the quite tractable expression 

5X ~-XTy"(C+XM)x b /[y b x b + Xy"Mx b l (5.26) 

or equivalently, after multiplying throughout by A and using (5.19 a ), 

SX ~ -A 2 r- y H a (C + XM)x b / yf (AW - K)x b . (5-27) 

This can be computed easily for any given structure, so giving both the direction and magnitude 
of the change to be expected in each closed-loop pole as a result of a specified delay. 
Examination of that pole which goes unstable most quickly therefore allows an estimate to be 
made of the critical delay r* at which the system will first be expected to become unstable. 


Note that Eqs. (5.26) and (5.27) are invariant under scaling of the left and right generalized 
eigenvectors. It is therefore not necessary to carry out any normalization procedure on these 
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vectors before computing the eigenvalue perturbation using this expression. The only 
computational steps required are therefore: 

(a) Construct the matrices M , M and C using (2.9), (3.5) and (3.6), respectively. 

(b) Solve for the eigenvectors, using (5.24) and (5.25), for the particular X of interest. 

(c) Compute the scalar expression (5.26) or (5.27). 

In most practical applications, the physical structure to be controlled will be very lightly damped. 
Furthermore, the AVA damping will also generally be chosen to be quite low, as will be further 
illustrated by the examples in the next section. It is therefore of interest to consider how delays 
affect the closed-loop eigenvalues in this case. Many of the terms in (5.27) then simplify 
considerably: in particular, we have that X~j(d, y a = x b is approximately real, and M«C. 

We can therefore rewrite (5.27) as 

5X - -Q) 2 T-\ T b C\ b / x T b (co 2 M + K)x b . (5.28) 

All terms in this approximation are real, so the eigenvalue perturbation for lightly-damped 
structures is too. Thus, feedback delays mainly affect only the real parts of the eigenvalues of 
such systems. 

It is interesting to note that feedback delays will not destabilize all modes of the composite 
system. Instead, what typically occurs is that the low-frequency modes are stabilized, i.e. shifted 
farther into the left half-plane, while it is only the high-frequency modes that are destabilized. 
The crossover between these two regimes is just, for the case of a single degree-of-freedom 
AVA, the natural frequency of the AVA virtual system. (The generalization to multi-degree-of- 
freedom AVAs is quite straightforward.) This can be shown by means of (5.28) as follows: the 
matrix M is positive definite, so co 2 M + K certainly is also. Thus, the denominator of this 
equation is guaranteed to be positive. The sign of 5X will therefore just be the opposite of the 
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sign of the term x T b Cx b . But, from (2.6) and (5.6), this can be written in terms of the 
displacements of the physical and virtual structures as 

x T b Cx b =-q T B® v diag(col)T\ v . ( 5 - 29 ) 

Now, the input to the AVA is just, by (2.3), u v = B T q. Regarding this as a lightly-damped single 
degree-of-freedom harmonic oscillator, its output r] v will be approximately in phase with its 
input at frequencies below co v , and 180° out-of-phase with it at frequencies above it. But u v is 
itself 180° out-of-phase with the term S r q , as a result of the double derivative. Thus, (5.29) will 
be positive at low frequencies and negative at high. This demonstrates that delays do indeed 
stabilize low-frequency modes while destabilizing the higher-frequency ones. 


It should be noted that the results presented here are based on an assumption of distinct closed- 
loop eigenvalues. If the composite system has repeated frequencies, the first-order eigenvalue 
perturbation expression (5.18) no longer applies as it stands. Fortunately though, similar results 
can be still obtained by means of the methods of [20]. The main complication in this case arises 
from the fact that the eigenvectors corresponding to a set of repeated eigenvalues are not unique; 
rather, they span a subspace of permissible eigenvectors. There is therefore a degree of non- 
uniqueness in the definition of the eigenvectors to be used in the first-order eigenvalue 
perturbation equation. However, it is still possible to define a set of differentiable eigenvectors 
[20] corresponding to the repeated eigenvalues, and these unique vectors can be used to obtain 
well-defined eigenvector perturbations. Thus, the results obtained here generalize directly to the 
repeated-eigenvalue case. 

A particularly simple, although somewhat conservative, lower bound on the critical delay can be 
obtained by considering the definiteness of the composite mass and stiffness matrices. As 
already noted, these matrices are not symmetric once a delay is introduced into the feedback 
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loop. However, their eigenvalues can be shown to be, to first order, equal to those of their 
symmetric parts 

M s = \{M X + M t x ) = M + j(M + M t ) 

and (5.30) 

C, = i(C, + C t r ) = C+|(C+C r ). 

A sufficient, but not necessary, condition for closed-loop stability is that these matrices both be 
positive (semi-)definite. Thus, a lower bound on the critical delay is the delay for which either 
M, or C s first has an eigenvalue which crosses the imaginary axis. For lightly-damped systems, 

it should be noted that the eigenvalues of the nominal damping matrix C will be considerably 
smaller than those of the nominal mass matrix M . Furthermore, the damping perturbation C is, 
by (5.5) and (5.6), significantly smaller than the mass perturbation M. The damping matrix will 
therefore become indefinite for a far smaller delay than will the mass matrix; it is consequently 
the only matrix that need actually be tested. 


From (2.9) and (5.6), the symmetric pseudo-damping matrix C s is given as 


C,= 


(C 0 

0 diag(2£ Vi co vt )) 


^ T f 


0 B<t> v diag(0)l)^ 

diag(o) 2 vl )<t>l B r 0 


(5.31) 


Now, this matrix becomes indefinite when one of its eigenvalues equals zero, i.e. when 
det{CJ = 0. Making use of the same type of determinantal identity [12] that led to (4.4), this can 

be rewritten as 


det{C - T 2 • B& v diag(-^y-)<t>l B t } = 0. 


(5.32) 


This matrix can be regarded as the nominal damping matrix C perturbed by a term proportional 
to the small quantity r 2 . Standard eigenvalue perturbation results [11] [21] can therefore be used 
to predict, to first order, the value of T for which one of its eigenvalues will equal zero. Suppose 
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C has eigenvalue /i and corresponding eigenvector v, normalized so as to give ||v|| 2 = 1. (Note 
that the symmetry of C implies that v is both the left and right eigenvector corresponding to fi .) 
Then, a first-order approximation to the perturbation in this eigenvalue is just 




n 


til 



(5.33) 


Thus, the delay at which this eigenvalue becomes zero is approximately given as 


r. - 


i 




cot 


(5.34) 


v T B® v diag(— r~)& v B v 

°Svi 


Considering all eigenvalues of C gives n expressions of this form; the smallest of these is the 
estimated delay at which indefinite damping is obtained. This expression provides insight into 
the way in which the AVA gain-like parameter <b v affects delay robustness: as ||d> v || increases, so 

does the denominator of (5.35). The delay at which C ceases to be positive definite therefore 
decreases, i.e. the closed-loop system becomes less tolerant to feedback delays. This is certainly 
physically reasonable. 


Equation (5.34) requires only the solution of an (n x n) symmetric eigenproblem, and so is quite 
attractive computationally. It becomes even simpler if the original physical structure is given in 
modal form. We then have that C = diagil^coj, so its eigenvalues are just {2 and its 
eigenvectors are of the form ((),••■, 0,1,0, •••,0) r . The critical delay (5.34) therefore reduces to 


T ** 4 • 

S J 


1 


C co 


p T diag(^ L )p 

bvr 


(5.35) 


where p = and is the f 3 * row of B. 
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6. Examplpft 

Tl.e .suits derived in the previous sections win now be iUust rated by app.ication to flexible 

structures of progressively increasing complexity. The first system considered is a single degree- 

of-freedom one, i.e. a physical structure consisting of a single mass, spring and dashpot. In this 

appl.cat.on the disturbance source, prima^ output sensor and AVA sensor/ac.uator pair are all 

necessarily pos.Uoned on the single mass, and so are cettainly coincident. A modal mode, of a 

cantilever beam is then studied, in order to investigate the parties of an AVA when applied to 

a system w„h multiple degrees-of-freedom. The AVA sensor/actuator pair is first taken to be 

coincident with either the single disturbance source or primary output sensor; multiple 

disturbances are then considered; and then the non-coinciden, case is investigated. Finally, the 

effects of feedback delays are studied for both the mass-spring-dashpo, and cantilever beam 
systems. 


Sing[g ^ Pggrgg~of~Freecf^rr| System 

Consider a single degree-of-freedom physical structure with a natural frequency of V2 md/s and a 

P ng ratio of 0.01. An AVA is required to introduce a transmission zero into this system at 

dte disturbance frequency of 0.6 rad/s; as the system is coincident, this is achieved simply by 
selecting m. = 0.6 rad/s as the natuml frequency „f the virtual system. 


The role of the virtual damping parameter f. can be observed in * closed-loop frequency 
response plots (Fig. „ that are obta^ed for f, =0.001,0.1 and 0, (and *. = 1 . 5 ); ,he „pe„., 0 op 
plot ,s also g.ven for comparison. A choice of f. = 0.1 thus gives a good compromise between a 
transmission zero bough that is deep and one tha, is wide, so providing considerable vibration 

::„d t of frequencies - ume ~ ^ - — ** - dam Pi „g 

er armomc excitation at the disturbance frequency indicate that f, = 0. 1 will 

also give fas, transient decay and small steady-state ereor; this damping rado is therefore chosen 
as the design value. 
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The only parameter remaining to be selected is the scalar <D V ; this will be chosen in this example 

so as to minimize the condition numbers of the closed-loop poles, using (3.8) and (3.9). Fig. 3 
shows the corresponding plots; a choice of <h v =0.8 yields closed-loop poles which are as robust 

as possible, and completes the AVA design procedure. 

6.2 Coincident Cantilever Beam 

As an example of a structure with multiple degrees-of-freedom, consider the transverse 
vibrational dynamics of a uniform cantilever beam of length 25 m, width 0.1 m and depth 0.01 
m, with material properties those of aluminum (density 2.7x10^ kg/m^, Young's modulus 
7.0xl0l0 N/m^). The primary displacement sensor is located at the free tip of the beam, and a 
harmonic disturbance source operating at 0.4 rad/s is positioned 15 m from the built-in end. A 6- 
mode model of this structure has natural frequencies between 0.0827 and 7.0213 rad/s; this 
model will be studied here, with a damping ratio of 0.005 assumed for all modes. 

An AVA is to be designed to minimize the effects of the disturbance signal on the measured 
output. A virtual damping ratio of = 0.10 will be selected in all the configurations studied, 

based on the observations of the last sub-section. As a first case, Fig. 4 shows the frequency 

responses that are obtained by placing the AVA at the disturbance source; the three plots are for 
values of the gain-like quantity G> v of 0 (open-loop), 5 and 10. It can be seen that increasing the 

gain increases both the depth and width of the virtual zero trough, as should be expected. 
Furthermore, the three higher-frequency structural zeros clearly do not shift from their open-loop 
values; this agrees with the result proved earlier for coincident systems. Fig. 5 then shows the 
variation of the LQR-like objective function J of (3.7) with <!>„; a control weighting of p =0. 15 
was used. A choice of <I> V = 3.3 is thus optimal in this sense, and can be taken as the design 
value. 
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Fig. 6 is analogous to Fig. 4, but with the AVA now coincident with the primary output sensor 
rather than the disturbance source. The basic properties of these two configurations are quite 
similar, although the details can be seen to differ somewhat. Fig. 7 then shows the closed-loop 
output obtained in the latter case in response to the harmonic disturbance excitation, which has 
unit amplitude; a fairly high value of <I> V = 7 is used in order to ensure good output regulation. 

The effectiveness of the AVA in absorbing the sinusoidal excitation can clearly be seen. 

Finally, the use of a single AVA to absorb two distinct disturbance signals will be illustrated. A 
second harmonic disturbance acting at 2.0 rad/s is applied to the beam at a location 7.5 m from 
the built-in end, and a two degree-of-freedom AVA positioned coincident with the primary 
output sensor. The natural frequencies of the virtual structure in this case are just the two 
disturbance frequencies, and the virtual modal parameter O v is now a (1 x 2) vector. Fig. 8 

shows the closed-loop frequency responses from each of the two disturbance stations to the 
measured output that result for a choice of 4> v = (7 5). Fig. 9 then gives the corresponding 

closed-loop time response obtained if the first disturbance has a peak amplitude of 1 and the 
second 2. By comparison with the open-loop response (which clearly exhibits the effect of the 
high-frequency input: see Fig. 7), the AVA can be seen to again produce very acceptable 
rejection in this multi-disturbance case. 

6.3 Non-Coincident Cantilever Beam 

We shall now briefly illustrate of the design of a single-mode AVA with a sensor/actuator pair 
which is not coincident with either the disturbance or output. Consider the same cantilever beam 
as before, with the primary sensor still at the free tip. However, the AVA is now located 10 m 
from the built-in end, and a disturbance source of unspecified frequency is positioned at 15 m. 
Equation (4.9) is then used to determine the required AVA natural frequency and damping ratio 

as the disturbance frequency ranges over all possible values; in all cases, a damping ratio of 0. 1 
is specified for the virtual zero. Fig. 10 shows the resulting virtual frequency £O v versus 
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disturbance frequency for a choice of <!>, = 7, and Fig. 1 1 gives the corresponding £ v plot. It can 
be seen that the only large perturbations occur in the vicinity of the structural zero near 1 rad/s, 
as predicted. Finally, Fig. 12 shows the resulting closed-loop frequency response plot produced 
by this non-coincident controller, together with the open-loop plot for comparison. The AVA 
can be seen to have indeed introduced a zero at the desired frequency. 


6.4 Effects of Delays 

The remaining six plots illustrate the closed-loop stability properties of the two structures 
discussed above when a small delay is present in the AVA feedback loop. The mass-spring- 
dashpot system AVA is chosen in this case to have virtual modal parameters co v =1.0 rad/s, 
£ v =0.02, and either = 1 (referred to as the low gain case) or <I> V =3 ( high gain). The 
resulting loci of closed-loop poles as a function of the delay r are given in Figures 13 and 14. In 
both plots, the poles obtained for x = 0 are indicated by 'o's. It can be seen that the lower- 
frequency AVA mode is stabilized in both cases by the introduction of the feedback delay, while 
the higher-frequency mode is destabilized. This behavior is indeed of the type predicted by 
(5.29). Furthermore, the low gain system becomes unstable for a delay of 0.076 s and the high 
gain one for 0.048 s. This agrees with the expectation that the critical delay will decrease with 
increasing AVA control gain. Figure 15 provides further information on this point by plotting 
against control gain <I> V the actual critical delay (solid curve); that predicted by the eigenvalue 

perturbation expression (5.27) (dashed plot); and the simple conservative bound obtained from 
the damping matrix (5.35) (dotted plot). The very close agreement between the true critical 
delay and that predicted by (5.27) can be easily seen, as can the general trend of decreasing delay 
with increasing <!> v . It can also be seen that the values of critical delay are considerably smaller 

than the periods of the two modes of this system; this confirms that the first-order assumption 
made in (5.1) is indeed acceptable. 
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Finally, Figures 16 - 18 are analogous to Figures 13 - 15, but for the 6-mode cantilever beam 
model discussed above. The AVA in this case is positioned at the free end of the beam, with 
modal parameters taken to be exactly as for the mass-spring-dashpot delay studies. It can be 
seen from Figures 16 and 17 that the low-frequency modes are again stabilized by delays and the 
high-frequency ones destabilized, as predicted. In addition, it is important to note that the 
highest mode is noi actually the critical one for this system. It is therefore necessary to test 
(5.27) for each of the modes with frequencies above the AVA frequency. (If a multi-degree-of- 
freedom AVA were employed, modes with frequencies lying between the AVA frequencies 
would also have to be checked.) Lastly, Figure 18 demonstrates the same close agreement 
between true critical delay and that predicted from the eigenvalue perturbation results as was 
obtained for the mass-spring-dashpot system. The only significant difference is the slight 
discontinuity in the beam critical delay curve for a gain of approximately 4.3. The explanation of 
this is simply that a different mode becomes the critical one for gains higher than this value. 
This again reflects the fact that all higher-frequency modes are potentially destabilized by delays, 
and so should be tested by means of the perturbation expression (5.27). 
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7. Conclusions 

This paper has analyzed the properties of virtual passive controllers (active vibration absorbers), 
and of the closed-loop systems achievable by means of such compensation. In particular, it was 
shown that these controllers introduce additional transmission zeros into the closed-loop system; 
they therefore provide an extra degree of design freedom over that provided by state feedback, 
which cannot alter the zeros in any way. 

The control design problem of providing good vibration at a specified frequency, for instance 
that of a disturbance, is equivalent to requiring that a zero be added at this frequency. 
Techniques were developed for designing a virtual passive controller to achieve this goal; these 
methods were shown to be particularly straightforward in the case where the controller 
sensor/actuator pairs are coincident with either the disturbance sources or the primary output 
sensors. Methods were then derived for exploiting the remaining controller design freedom to 
achieve desired secondary objectives, for instance making the closed-loop poles as robust as 
possible. 

Finally, the effects on closed-loop stability of small delays in the feedback loop were 
investigated. It was shown that these delays tend to destabilize the higher-frequency modes 
while actually making the lower-frequency ones more heavily damped. Eigenvalue perturbation 
results were then derived which allow accurate predictions to be made of the critical delay at 
which instability first occurs. These points were illustrated by simple examples. 
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Figure 1 1-DOF system frequency responses 
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Figure 8 Frequency responses, two disturbances 







Figure 9 Time responses, two disturbances 















Figure 12 Frequency responses, non-coincident 
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Figure 13 1-DOF root locus 
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Figure 14 1-DOF root locus against delay, high gain. 
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Figure 16 Beam root locus against delay, low gain. 







m 

o 


m 

o 






o 








